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Abstract 

A new and simple statistical approach is performed to calculate the parton distri- 
bution functions (PDFs) of the nucleon in terms of light-front kinematic variables. 
We do not put in any extra arbitrary parameter or corrected term by hand, which 
guarantees the stringency of our approach. Analytic expressions of the x-dependent 
PDFs are obtained in the whole x region [0,1], and some features, especially the 
\cw-x rise, are more agreeable with experimental data than those in some previous 
instant-form statistical models in the infinite-momentum frame (IMF). Discussions 
on heavy-flavored PDFs are also presented. 
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1 Introduction 



One of the goals in fundamental physics is to search for the detail informa- 
tion of the nucleoli structure. The parton model, suggested by Feynman [I], 
and then immediately applied to the deep inelastic lepton-nucleon scattering 
process by Bjorken [2], proposes that the nucleoli is composed of a number of 
point-like constituents, named "partons" , which were afterward recognized as 
quarks and gluons. In the impulse approximation, the deep inelastic lepton- 
nucleon scattering can be viewed as a sum of elastic lepton-parton scattering, 
in which the incident lepton is scattered off a parton instantaneously and in- 
coherently. This is in accord with one property of QCD - asymptotic freedom. 
On the other hand, due to another property of QCD - color-confinement, the 
constituents of nucleon - quarks and gluons, have never been seen individually. 
The nucleon structure functions, in terms of the parton distribution functions 
(PDFs), are badly desired in hadronic study. However, due to the complicated 
non-perturbative effect, we still have difficulty to calculate them absolutely 
from the first principal theory of QCD at present. 

Various models according to the spirit of QCD have been brought forward, 
therein statistical ones, providing intuitive appeal and physical simplicity, have 
made amazing success $MflfH7ltMnMttMM%nHnKnnMtMlffl\ . Ac- 
tually, as can be speculated, with partons bound in the wee volume of the 
nucleon, not only the dynamic, but also the statistical properties, for exam- 
ple, the Pauli exclusion principle, should have important effect on the PDFs. 
Cleymans and Thews [3|H] , as pioneers, started with the transition rate of scat- 
tering in the framework of temperature dependent field theory and explored 
a statistical way to generate compatible PDFs. Mac and Ugaz [5] incorpo- 
rated first order QCD corrections (however, the perturbative term turned out 
to be a sizable fraction of the statistical term), and afterwards Bhalerao et 
al. [TT1ll2lll3|fn] introduced finite-size correction and got more fitting results; 
they both referred to the infinite-momentum frame (IMF). Bickerstaff and 
Londergan [6] interpreted the finite-temperature property as to mimic some of 
the volume- dependent effects due to confinement, furthermore, they discussed 
the theoretical validity of the ideal gas assumption in detail. Devanathan et 
al. [SfHyTU] proposed a thermodynamical bag model which evolves as a func- 
tion of x, and the structure functions they got are practicable for x > 0.2 and 
have correct asymptotic behavior as x — > 1; in addition, they parametrized on 
T and exhibited the scaling behavior. Bourrely, Soffer and Buccella [T5fl6lll7] 




developed a new form of statistical parametrization, allowing x-dependent 
chemical potential, and by incorporating QCD evolution they got indeed re- 
markable PDFs. Otherwise, Zhang et al. [T8lll9f20j constructed a model using 
the principle of detailed balance and balance without any free parameter, and 
the Gottfried sum they got is in surprisingly agreement with experiments. 
Alberg and Henley [21] tracked the detailed balance model for a hadron com- 
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posed of quark and gluon Fock states and obtained parton distributions for 
the proton as well as the pion. 



When dealing with physics in high energy region, light-front dynamics is a suit- 
able language [22] . It is known [23] that the impulse approximation fails when 
using instant-form dynamics in the nucleon rest frame, but works well when 
using it in the IMF or using light-front dynamics in an ordinary frame. How- 
ever, statistics in light-front formalism encounters some difficulties, therein 
the most fatal one is the unclarity about what the light-front temperature 
is and how it relates with the usual one in instant-form formalism (see, e.g., 
Refs. [24,25 26 27]). Consequently, the way of generalization from instant-form 
statistics to light-front statistics is quite speculative. Actually, even the gener- 
alization of the thermal dynamics and statistical theory from the system rest 
frame to a moving frame in instant-form formalism is also not so understood, 
and discussion on it has continued for a long period (For theoretical discussion, 
see Refs. [28129 30 31 32 TI331I34"] and references therein, and for recent numerical 
experiments, see, e.g., Refs. [351136] ). 




In order to avoid these tough problems, we start with instant-form statistical 
expressions in the nucleon rest frame, then perform transformation on them in 
terms of light-front kinematic variables. The analytic expressions of the PDFs 
we get are something different from those attained in other statistical models 
performing in the IMF [5llllfl2M13|[T4"] . The largest distinction is that, when 
x — > 0, the distributions of our light-flavored quarks (anti-quarks) do not tend 
to zero as theirs, but give a rise instead, which agrees with the experimental 
data better. 

Worthy to note that, our intention is only to illustrate whether the statistical 
effect is important and to which aspects of the nucleon structure it is impor- 
tant, not how well it matches the experimental results, so we do not make 
any effort to fit the experimental data intentionally. There is no arbitrary pa- 
rameter put by hand in our model, and all parameters are basic statistical 
quantities. Some of other statistical models can fit the experimental data bet- 
ter by introducing many free parameters, however, it weakens the stringency 
at a cost. In addition, our results naturally cover the whole x region [0,1], and 
the features of PDFs and structure functions at the boundary are of great 
interest in both theoretical and experimental study. 

The paper is organized as follows. In section 2, a brief description about the 
approach used in this paper is introduced, and the analytic expressions of 
the PDFs are presented. In section 3, numerical results and comparisons with 
experiments and other theories are illustrated. In section 4, the mass effect of 
the partons and the features of heavy-flavored PDFs are discussed. The last 
section is a short summary. 
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2 The statistical approach 



We assume that the nucleoli is a thermal system in equilibrium, made up of 
free partons (quarks, anti-quarks and gluons). In the nucleon rest frame, the 
mean number of the parton (denoted by /) is 

N f = J f(k°) d 3 fc , (1) 

where f(k°) satisfies the Fermi-Dirac distribution or the Bose-Einstein distri- 
bution 

with the upper sign for Fermion (quark, anti-quark), and nether sign for Boson 
(gluon); gf is the degree of color-spin degeneracy, which is 6 for quark (anti- 
quark) and 16 for gluon; /i/ is its chemical potential, and for anti-quark fig = 
—fi q , for gluon fig = 0. 

Instead of boosting above expressions to the IMF [5flllll2lll3|[T4"] , we transform 
them in terms of light-front kinematic variables in the nucleon rest frame. 
Before doing this, note that when doing the integration in Eq. ([!]), the on- 
shell condition k° = ^/k 2 + m 2 is needed, where k°, k = (fc 1 , A; 2 , fc 3 ) , m/ are 
the energy, 3-momentum and mass of the parton, respectively. Eq. (CTJ can be 
explicitly reexpressed as 

N f = J f(k°)5 (k° - ^(k^ + kl + mj) dk°dk 3 d 2 k ± . (3) 

We introduce the light-front 4-momentum of the parton k = (k + , k~ , kj_), 
where k + = k° + k 3 , k~ = k° — k 3 , h± = (k 1 , k 2 ), and k + = P + x = Mx, where 
x is the light-front momentum fraction of the nucleon carried by the parton 
and M is the mass of the nucleon. Hereby, the 5-function and the integral in 
Eq. © turn to 



5 (k° - ^(k 3 ) 2 + k 2 ± + m) 



2k°9{k°)S(k 2 -m) 



1 + 



k^ + m 2 



(Mx) 



9(x)5\ k~ 



k 2 ± + m 2 ' 

Mx 



and 



dk°dk 3 d 2 k { = -Mdk~dxd 2 k , . 



(4) 



(5) 



Then Eq. ((2]) becomes 



N f = J f{x, k ± ) d 2 k L dx 



(6) 
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where f(x, kjj is 



f(x, kj 



g f MV 



2(2tt) 



Mx+- 



TTT7 



1 + 



(Ms) 2 



0(z) 



(7) 



± 1 



On the trivial assumption that is transversely isotropic, we can integrate 
on |kj_| analytically and get 



/(*) 



g f MTV 
8tt 2 



Mx 



-2TLi 2 



m, 



Mx 



In 



l±e 



Mx+ 



2 
/ 

TTTr 



V 



/J 



as is mentioned above, the upper sign for Fermion and the nether sign for Bo- 
son. Li 2 (z) is the polylogarithm function, defined as Li 2 (z) = J2h=i z k /k 2 - Note 
that the expressions of the PDFs (Eq. (jSJ)) are different from those attained 
in the previous statistical models [3lllll5ir7ffTT1ll2|13|14j . 

In practice, the PDFs in a certain system should be constrained with some 
conditions. For example, in the proton, they are 



Uy 



[u(x) — u(x)] dx 



(9) 



dy = \d(x) — d(x)] dx = 1 , 
Jo 

[u(x) + u(x) + d(x) + d(x) + g(x)]x dx = 1 . 



(10) 
(11) 



Thus for the proton, we have three equations (191) (TTOT ) ffTTT) . and four unknown 
parameters T, V, fi u , [id (The mass of the proton, M = 938.27 MeV, is taken 
as given). So for a given T, the rest parameters V, /i u , [id can be determined 
by solving the equations. 



3 Results 



We perform our calculation for the proton, therefore the following results, if 
not specially stated, are all for proton. However, the method is absolutely 
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applicable to the neutron. For convenience, we just consider the u, d flavor 
and gluon, and take rrif = 0, which will be showed, in section 4, as a good 
approximation. 

In practice, we adopt a certain value of T at first, then numerically solve the 
equations to get V, fi u , fid- Subsequently, the PDFs can be obtained according 
to Eq. (jSJ), as well as the Gottfried sum 

^ = ri F|(x)-F 2 "(x) ^ = 1 2 ^ r _ g i ^ 
Jo x 6 6 Jo 1 J 

We find that at T = 47 MeV, the Gottfried sum Sq = 0.236, which agrees 
well with the experimental result Sq = 0.235 ± 0.026 [37]. We conclude that 
the temperature of proton is around 47 MeV, and V$ ~ 1.2 x 10~ 5 MeV~ 3 , 
fi u ~ 64 MeV, fid ~ 36 MeV. The following results are all given at To = 47 
MeV. 

Taking proton as a perfect sphere, we can calculate its radius ro from the 
volume Vq. We get r = 2.8 fm, which seems a little bigger than the practical 
value, possibly due to the oversimplified assumption of the uniform distribu- 
tion of partons. Worthy to mention that, the r we get, together with T , is 
close to what Mac and Ugaz [5] got with the consideration of first-order QCD 
corrections, and their fitted values are T = 49 MeV, r = 2.6 fm. 

The PDFs f(x) and xf(x) are illustrated in Fig. [Hand Fig. El respectively. 
In contrast to q(x) and q(x) — > as x — > in the previous statistical models 
without extra corrected term [3|4][5|7] . our trend makes a good improvement. 

In our model, the flavor asymmetry of the nucleon sea is naturally generated. 
d{x) — u(x) and d{x) / u(x) are shown in Fig. [3] and Fig. [H the former agrees 
well with the results of experiments and CTEQ parametrization [38] while the 
latter seems not. The model suggested by Zhang et al. [T8ll9ll20j also gives 
reasonable asymmetry of u and d without introducing any parameter, which 
is further discussed by Alberg and Henley [21]. And the feature of d(x) —u{x) 
in Refs. [T9l[2T] and d(x)/u(x) in Ref. [21] have similar behavior as ours. 

Furthermore, d(x)/u(x) is also shown in Fig. [5) Note that when x — > 1, 
d{x)/u{x) — > 0.55, which is different from the result of CTEQ parametriza- 
tion [38], but close to the prediction of the naive SU(6) quark model. 

The nucleon structure function F 2 (x) = 2xFi(x) = xJ2f e }f( x ), where ej is 
the charge of the parton of flavor /, is shown in Fig. [6J With the p-n isospin 
symmetry, i.e. u n (x, kj_) = d p (x, kj_), d n (x, kj_) = u p (x, kj_), u n (x, kj_) = 
d p (x, kj_), d n (x, kj_) = u p (x, kj_), g n (x, kj_) = g p (x, kj_), we can also obtain 
the structure function of the neutron F£(x). Ff (x) — F^x) and F^x) / F^x) 
are shown in Fig. Hand Fig. El We can see that F^ix) — F£{x) is quite agree- 
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Fig. 1. The calculated f(x) in our statistical approach. 



1.0 




Fig. 2. The calculated xf(x) in our statistical approach. 

able, while F^x) / F^x) is not. One interesting feature is the behavior of 
F£(x) I F^x) when x — > 1. In the naive SU(6) quark model it tends to 2/3, 
while in the SU(6) quark-diquark model it is 1/4 and in experimental observa- 
tion it is smaller than 1/2. Here our result seems close to the prediction of the 
naive SU(6) quark model again, however, it does not agree with the reality. 

From the above results, we find that our statistical method can successfully 
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Fig. 3. Comparison of our result with CTEQ result at Q 2 = 1 GeV 2 [38], 
E866/NuSea result at Q 2 = 54 GeV 2 [39] and HERMES result at < Q 2 >= 2.3 
GeV 2 |S] for d{x)-u{x). 
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Fig. 4. Comparison of our result with CTEQ result at Q 2 = 1 GeV 2 [38 
E866/NuSea result at Q 2 = 54 GeV 2 [39] and NA51 result 00] for d(x)/u(x). 
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Fig. 6. Comparison of our result with CTEQ result at Q 2 = 1 GeV 2 [38] for F 2 (x). 

describe the behavior of the "subtracted-terms" , such as d(x) — u(x) and 
F$(x) — -F^(x), but the "divided-terms", such as d(x)/u(x), d(x)/u(x) and 
i 7 f(x)/F 2 "'(x), can only match the trend of experimental results approximately, 
and the departure is especially large in the high-x region, where the valence 
parts of the PDFs dominate. This feature probably implies that an additive 
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Fig. 7. Comparison of our result with CTEQ result at Q 2 = 1 GeV 2 [38] for 
F£(x)-F?(x). 
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Fig. 8. Comparison of our result with CTEQ result at Q 2 = 1 GeV 2 [38] for 

statistics-irrelevant corrected term to the PDFs works, whereas more free pa- 
rameters and uncertainty should be introduced. Bhalerao et al. successfully 
reproduced most features of the PDFs and structure functions by introducing 
two extra corrected terms [TT|12|13fT3] , at the cost of two more free parameters 
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Fig. 9. d{x) — u(x) and d(x)/u(x) at different masses (unit: MeV) of u, d quark. 



and violating the p-n isospin symmetry. 
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4 Further discussions 



We have ignored the masses of the quarks and anti-quarks for simplicity. Nev- 
ertheless, mass effect can be taken into account without difficulty. Actually we 
have performed this calculation and found, as can be speculated, the correc- 
tion of mass effect to the light-flavored PDFs is very small. d(x) — u(x) and 
d(x)/u(x) with different masses of u, d quark are illustrated in Fig. [9j 

However, the mass difference between u and d quarks can generate the mass 
split between the proton and the neutron. The invariant mass square of the 
system is given by 

/ m 2 + k 2 ^ 



x 



m2 = lly—r^) ■ 

In the continuous condition, it is 



M 



E 

/ 



m. 



x 



f(x,k ± )d 2 k A 



dx . 



(14) 



Using the p-n isospin symmetry, we get 

1 



Mt - Mi 



x 



u p (x) + u p (x) - d p (x) - d p (x) 



dx . 



(15) 



In PDG 2006, m u = 1.5 ~ 4 MeV, m d = 3 ~ 7 MeV, A pn = M n - M p = 
939.565 - 938.272 = 1.293 MeV. 

When we use the mean value m u = 2.25 MeV, rrid = 5 MeV, we get A pn = 
0.664 MeV from Eq. ([TBI) , and when we use the extreme value m u = 1.5 MeV, 
m<2 = 7 MeV, then A pn = 1.557 MeV. The result seems rather agreeable and 
it reinforces the reasonableness of our approach. 

We have only calculated the light-flavored PDFs, however, the heavy flavors, 
such as s, c, b, t quarks and the corresponding anti-quarks, can be treated 
in the same way. Since, in nucleon, the valence numbers of them are zero, 
the chemical potentials of them must all be zero. It leads directly to three 
following conclusions: 

Firstly, except for their masses, the heavy-flavored PDFs have no additional 
free parameters than the light-flavored ones. That is, if their masses are used 
as inputs, their PDFs can be uniquely determined by the parameters T and 
V, which have already been fixed in the previous light-flavored condition. So 
the difference between heavy-flavored PDFs only comes from the difference of 
their masses. 



12 



10 




d s(m s =100)-XN 



' T 0.2 f OA 
s(m s =200) S (m =300) 



0.6 0.8 



1.0 



Fig. 10. Comparison of s(x) at different mass (m s = 100,200,300 MeV) with the 
light-flavor PDFs at T and V . 

Secondly, the quark and anti-quark of the same heavy flavor have just the 
same distribution. For example, the condition s(x) = s(x) holds in the whole 
region x G [0, 1], so that the s, s asymmetry in the nucleon [12] does not come 
from the pure statistical effect. 

Thirdly, from Eq. ([8]), we can see that f(x) decreases when rrif increases. s(x) 
with different m s at To and Vq, together with the light-flavored PDFs, are 
illustrated in Fig. [TDJ and it shows that when m s < 100 MeV the contribution 
of s quark is considerable, and when m s > 200 MeV it is minor. Therefore, 
the contribution of heavier flavors is negligible. Calculation also indicates that 
s quark contributes less than 7% both to the total light-front momentum 
fraction x (see Eq. ( IrTj) ) and to the total invariant mass square of the system 
at T and V (see Eq. (ITijl ). 



5 Summary 



We preform a simple statistical approach and obtain analytic expressions of the 
parton distribution functions in terms of light-front kinematic variables in the 
whole x region [0,1]. The low-x behavior of these parton distribution functions 
is different from those in some previous instant-form statistical models in 
the infinite-momentum frame and our results are more close to the reality. 
There is no arbitrary parameter or extra corrected term put by hand in our 
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model, which guarantees the stringency of our conclusion. Several features of 
the parton distribution functions and structure functions of the nucleon are 
compared with the results of experiments and other theories. Calculation of 
the mass split between the proton and the neutron is also performed. We have 
further discussions on the influence of the heavy flavors. All of these show that 
although the statistical effect is not everything, it is very important to some 
aspects of the nucleon structure. 
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